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ABSTRACT 


This thesis explores the application of the Weibull 
distribution to the Cost Effectiveness Analysis 
Spreadsheet Model (CEAMOD) for the Aircraft Engine 
Component Improvement Program (CIP). The current model 
assumes the exponential distribution for all unscheduled 
failures. This thesis explores whether the Weibull 
distribution can be used in the CEAMOD. This process is 
pg2rformed through the use of a simulation program and 
spreadsheet analysis. The thesis first examines past 
concerns that the use of the Weibull will create cyclic 
annual numbers of failures of a component over an engine's 
life. It then examines the empirical distributions of the 
number of failures per year generated from simulations for 
a range of Weibull parameters and compares these with the 
Poisson distribution resulting from the current CEAMOD 
assumption of exponential times between failures. 
Finally, it considers how life cycle costs considered by 
the CEAMOD will change if the Weibull distribution is 
assumed. The conclusion is that the nature of the 
distribution of the annual number of failures should be 
studied further before the Weibull is incorporated into 


the CEAMOD. 
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I. INTRODUCTION 


A. BACKGROUND 

This thesis analyzes the input and methodology of the 
statistical distributions applied to the General Electric 
Aircraft Engines Cost Effectiveness Analysis Spreadsheet Model 
(GE CEAMOD). 

The GE CEAMOD is a spreadsheet program, originally 
developed by Pratt and Whitney, and later refined by General 
Electric for use in analyzing the cost effectiveness of 
aircraft engine change proposals [Ref. 1], The current 
release. Version 2.1, is written for use with the Microsoft 
EXCEL 4.0 program. Both the Navy and the Air Force use the 
CEAMOD program in conjunction with their prime contractors to 
determine whether an Engineering Change Proposal, as submitted 
by the manufacturer, is cost effective over the proposed life 
of the engine/component. 

The CEAMOD program is used in support of the aircraft 
engine Component Improvement Program (CIP). The CIP has the 
following three objectives [Ref. 1]: 

1. To maintain an engine design which allows the maximum 
aircraft availability at the lowest cost to the 
government. 

2. To correct, as rapidly as possible, any design 
inadequacy which adversely affects safety of flight. 

3. To correct any design inadequacy which causes 
unsatisfactory engine operation or adversely affects 
maintainability and logistic support in service. 
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The CEAMOD program supports the CIP by allowing engineers 
and program managers to analyze the expected life cycle costs 
differential between a current component and a proposed change 
in that component. Because the life of a component is 
typically many years, all future expected costs or savings are 
expressed in net present value terms. Since many of the 
factors included in the analysis are estimates, the final 
program output is subject to interpretation. In addition, 
there are other "softer" factors that are considered as well. 
These include: safety, increased mission effectiveness, and 
decreased overall downtime for maintenance and repairs. The 
CEAMOD does not analyze these factors, but none-the-less, they 
are critical to the overall decision of whether to incorporate 
the modification. 

One factor that has been under discussion with the 
program developers and users is the assumption in CEAMOD of a 
constant component failure rate. By assuming that the failure 
rate of a component is constant, the time between failures can 
be described by the exponential probability distribution. 
There is a belief that assuming this distribution ensures the 
CEAMOD results, which show fairly constant numbers of failures 
per year, are representative of reality. However, the 
designers and users of CEAMOD have proposed that the Weibull 
distribution may reflect the reality of the times between 
component failures more accurately. The concern is that due 
to the Weibull's ability to reflect changing failure rates. 
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the CEAMOD program may become too complex [Ref. 9]. In 
addition, the question remaining to be investigated is: "can 
the Weibull be successfully introduced into the CEAMOD program 
without creating a cyclic nature to the number of component 
failures over time?" 

B. OBJECTIVE 

The objective of this thesis is to analyze one of the 
specific assumptions of the CEAMOD program. This assumption 
has to do with the current and proposed unscheduled failure 
rates for the component in question. The current program 
assumes that the time between failures of a component is 
distributed according to the exponential probability 
distribution. If the time between failures behaves according 
to the Weibull distribution, how are the number of failures 
per year affected? There are two relevant issues with regards 
to the objective, they are: 

1. Can the number of failures per year be described by 
a simple probability distribution? 

2. Are there any significant differences in the overall 
life-cycle costs for a component when the 
exponential distribution is replaced by the Weibull? 


C. SCOPE 

The scope of this thesis is limited to an analysis of a 
single component's change and how the use of the Weibull 
distribution instead of the exponential affects the expected 
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life cycle costs. This thesis will only focus on the expected 
costs of unscheduled maintenance (i.e. unexpected failures) 
and not on the costs associated with planned maintenance. The 
trial data used was provided by General Electric and was the 
example previously analyzed in Rau's thesis [Ref. 4]. 


D. METHODOLOGY 

Due to the nature and limitations of the CEAMOD program, 
a simple substitution of the number of failures generated 
assuming the Weibull distribution was not possible [Ref. 9]. 
The exponential distribution provides for a constant failure 
rate. Thus, if a component reaches a certain point in its life 
where its failure rate is constant, it has no greater or 
lesser chance of failing in the next instant of time for any 
time during its life. The number of failures per year then has 
the Poisson probability distribution with the constant failure 
rate as its mean number per year. The Weibull distribution 
allows for a changing failure rate. The number of failures 
per year is no longer Poisson distributed. Thus, before the 
Weibull can be incorporated into the CEAMOD, the nature of the 
distribution of the number of failures per year must be 
understood. 

To study the distribution of the number of failures per 
year, assuming Weibull times between failures, a simulation 
program was written in the SIMAN simulation language. A 
version of the CEAMOD was also developed to accept the SIMAN 
data so that a financial analysis of the effects from the 
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Weibull distribution could be conducted. 


E. ORGANIZATION OF THE THESIS 

Chapter II provides a background of the current CEAMOD 
progratn and analyzes its methodology and assumptions. Chapter 
III is an in depth explanation of the proposed changes and 
their expected effects on the program. Chapter IV describes 
the development of the SIMAN program and the financial 
spreadsheet. Chapter V is an analysis of the model output. 
Chapter VI contains a summary of the research efforts and its 
conclusions. 
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II. BACKGROUND OF THE CEA MODEL 


A. DEVELOPMENT OF THE CURRENT MODEL 

The GE CEAMOD program is an EXCEL-based spreadsheet 
program currently used by the Navy and Air Force to estimate 
the future costs of aircraft engine component modifications. 
This program was originally developed as a mainframe program 
by Pratt and Whitney at the request of the Air Force. General 
Electric later converted it to a PC based Lotus 1-2-3 
spreadsheet and then to its current configuration as an EXCEL 
spreadsheet [Ref. l]. 

The Navy and Air Force have adopted the current model as 
the primary analysis tool in determining whether a 
modification to an engine or group of engines is cost 
effective from a life-cycle point of view. As noted in the 
Reeves thesis [Ref. 1], the model does have its imperfections 
and "illogical" steps, but it is currently the most effective 
model to use in determining long term cost effectiveness. The 
CEAMOD Users Group meets twice a year to discuss and approve 
improvements to the model. Version 2.2 was approved at the 
meeting held on 29 November 1994. 

B. ASSUMPTIONS OF THE CEAMOD PROGRAM 

There are a number of assumptions that have been 
incorporated into the current CEAMOD program. The primary 
assumption is that it will incorporate only one hardware 
change at a time. The significance of this assumption is that 
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a proposed component modification is considered to be totally 
independent of any other component. Thus, one modification 
has no bearing on any other. Although this may be true in 
Q 0 ]frain cases, it is not necessarily true in all cases. An 
aircraft engine is a complex assembly and the reliability of 
one component can and does have an impact on other components. 
This assumption was made primarily due to cost considerations 
in the development of the CEAMOD program. 

The costs usfed in the model to calculate the Net Present 
value (NPV) are primarily maintenance and operational costs. 
Also included in the model are the "kit" costs and labor 
required to make the change. When a component modification is 
proposed and the data is analyzed using the CEAMOD, the 
expected life-cycle costs associated with the current 
configuration are compared to those of the proposed changes. 
Figure 2.1 is a chart outlining which costs are included and 
which are not. 

Based upon these assumptions, the model is designed to 
ONLY include those logistics costs directly related to 
maintenance, operations and kit installation. All of the 
costs of developing the modification and program management 
are disregarded. For the most part it is a logical conclusion 
to not include the sunk costs of R&D and D&T, but if there is 
a significant life cycle cost differential between the current 
configuration and the proposed configuration for engineering 
data and program management these factors should be included. 
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COST CATEGORY 


COSTS INCLUDED IN MODEL 


Component Modification Cost 



Post-Modification Material 
Costs 


Cost of materials for 
maintenance AFTER modification 


Post-Modification Labor Costs Cost of labor for maintenance 

AFTER modification 


Cost of re-writing and 
implementing new documentation 


Cost of modifying or creating 
new support equipment 


Cost of purchasing new 
inventory 


Technical Documentation 
Modification costs 


Post-Modification Support 
Equipment Costs 


Post-Modification Inventory 
Costs 


Post-Modification Fuel Costs Cost of change in fuel usage 


Cost savings due to prevention 
of catastrophic aircraft loss 


Post-Modification Aircraft 
Loss Costs 


COSTS NOT INCLUDED IN MODEL 


Research and Development Costs Cost to develop modification 


Funds spent to ensure that the 


Design and Test Costs 


Engineering Data Costs 


Program Management Costs 


component performs as designed 


Cost to develop eng stds for 
new components 


Administrative costs related 


to the program management 

Table 2.1: CEAMOD cost Indusion/exclusion table 
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C. LIMITATIONS OF THE MODEL 

Although the CEAMOD program is quite complex, there are 
some inherent limitations. One critical limitation is the 
spreadsheet's ability to incorporate only one component change 
at a time. This limitation is primarily due to the program's 
ability to accept only one failure rate. This may be 
overcome, at least partially, by combining several component 
changes into one "complex" change. This process of combining 
failure rates is unrelated to the use of the program, and must 
be calculated offline [Ref. 1] . Another possible way to 
address the problem is to incorporate one component at a time 
in successive runs of the model. The iterative process might 
be to shift the modified configuration data to the "current" 
column and the next modification to the "proposed" column. By 
doing this for all components, the overall cost of all 
modifications can be estimated. This process of incorporation 
must be handled with care, because the order of incorporation 
may affect the outcome. Some orderings may increase overall 
costs because of interactions with other parts of the engine. 

Another limitation is related to the input data developed 
by the contract or/program manager. The data that is input 
into the program, both for costs and component reliability 
must be developed offline. Much of the data is estimated in 
the early stages of the program rather than being known 
precisely. Experience has shown that reality is different than 
the estimate. Thus, this program should not be used only at 
the inception of the program, but should be run at any time 
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there are significant changes to the relevant data. If the 
actual failure data of new components in service does not meet 
estimates, there may be a need to re-evaluate the modification 
policies. As a program develops, estimates of reliability 
and costs should be continually updated and entered into the 
program. Even small changes to reliability, cost or usage 
estimates may have a large effect on the total life cycle 
costs. These changes may lead the program manager to search 
for more reliable or cost-effective components. 

A significant limitation of this program is the 
assumption of the exponential probability distribution for 
times between failures. It will be the focus of this thesis. 
The assumption of the exponential distribution results in 
expected failure rates being constant [Ref. 4]. The 
consequences of assuming the exponential distribution for 
times between failure is that the number of failures over a 
specified period of time is Poisson distributed with a mean 
which is the product of the failure rate and the specified 


length of time. 

The input menu of the CEAMOD program requires a specific 
number to reflect the expected failure rate for both the 
current configuration and the proposed change. The rate is 
expressed as the expected number of failures per 1000 flight 
hours. The definition of "failure" is intentionally kept 
vague, and must be determined by the contractor/program 
manager. In the operating environment a failure may be defined 
as the time at which the component will no longer function, or 
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may be defined as the point at which it has failed to meet 
certain operating parameters. 

The assumption of a constant failure rate has been a 
topic of discussion for several years by the CEAMOD Users 
Group. The reason is that the reality of component failures 
is that they tend to decrease as the "burn-in" period passes, 
stabilize during the majority of the expected life, and 
increase in the later "wear-out" period. This phenomenon can 
be illustrated by the well-known "bathtub curve" [Ref. 10] . 
Figure 2.2 is an example of a failure rate bathtub curve. In 
those periods of decreasing and increasing failure rates, the 
time between failures are best estimated by Weibull 
distribution [Ref. 5]. How to incorporate the Weibull 
assumption into the CEA model has been a question which the 
users' group has puzzled over for several years. The purpose 
of this thesis is to attempt to determine how to do that 
incorporation. The specifics of the Weibull distribution 
will be discussed in Chapter III. 
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Figure 2.2; A failure rate bathtxib curve 
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III. PROPOSED CHANGES 


A. THE WEIBULL DISTRIBUTION 

The Weibull Distribution is named after a Swedish 
Physicist named Waloddi Weibull, who in 1939, developed it to 
represent the breaking strength of materials [Ref. 5]. Since 
then it has become widely used in many areas of analysis, 
particularly those related to the modeling of reliability of 
electrical and mechanical components [Ref. 5]. The Weibull 
distribution most frequently provides the best fit for the 
type of failure data experienced in the gas turbine industry 
[Ref. 5] . Thus, the consideration of using the Weibull 
distribution for times between failures in the CEA model has 
been an issue for the last several years. 

The probability density function (PDF), f(t,-11,13), of the 
Weibull distribution is represented by the following equation 
for values of t^O, ri>0, and |3>0 [Ref. 8] : 


t p 

n 


(3.11 


The specifics of the Weibull distribution are explained 
below, preceded by definitions of its parameters: 

1. BetaO) : Beta is commonly referred to as the "slope" 
parameter. The value of Beta O>0) determines which member of 
the family of Weibull curves that best fits the data being 
analyzed. The type of failure may be any of those represented 
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by the bathtub curve. "Slopes" or (3 values less than one 
correspond to decreasing failure rates, typical of those found 
in the infant mortality segment of the bathtub curve, 
"slopes" or (3 values equal to one represent "random" failure 
where the failure rate is constant, and "slopes" or p values 
greater than one correspond to increasing failure rates, 
typical of those found in the wear-out segment of the bathtub 
curve. A Weibull distribution with a beta of one is the same 
as the exponential distribution. 

2. Eta(ri): Eta defines the "characteristic life" of the 
Weibull distribution. The characteristic life of a component, 
as represented by eta, represents the value at which 63.2% of 
all failures will occur regardless of the value of beta. The 
value of l/h could be defined as the failure rate at the 
characteristic life in the Weibull distribution. The value of 
eta is the mean for the Weibull distribution (and, of course 
the exponential) when the value of beta is one, but is not the 
tnean when the value of beta is less than or greater than one. 
Figure 3.2 shows graphically the shape of the PDF for several 
values of beta. 

The hazard function describes the instantaneous failure 
rate. The hazard function for the Weibull distribution is 
represented by the following equation for value of t^O, h^O 
and p^O [Ref. 8] : 


]3(t;ri,P) = {^)^pt^''. 


(3.2) 
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Figure 3.2: Weibull PDF curves 


When (3 = 1 we have the exponential function, the hazard 

function becomes: 


Thus, because the failure rate equals l/n, and n equals the 
mean for an exponential distribution, the hazard rate is 
simply defined as one over the mean. 
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The mean of the Weibull distribution for values of beta 
less than and greater than one is given by [Ref.7]: 


P = 


ri 1 


(3.4) 


where r(.) is the Gamma Function. The variance of the Weibull 
distribution is given by [Ref. 7]: 


2 





(3.5) 


A comparative example showing the mean and variance for 
different values of beta is provided in table 3.1: 


p 

h 

t 

h(t) 

MTBF 

2 

.5 

50 

50 

.01 

100 

50000 

1 

50 

50 

.02 

50 

2500 

2 

50 

50 

.04 

44.3 

536.5 


Table 3.1: Comparison of means and variances 


Table 3.1 demonstrates that even though each example has 
the same t \ value (characteristic life) the actual mean values 
are very different. With a (3 = 1, or an exponential 
distribution, the mean is equal to the h value. But, when the 
beta value is less than one the mean is much higher and when 
the beta value is greater than one, the mean is less than the 
p value. The variance also decreases as the value of beta 
increases. 
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B. DETERMINATION OF WEIBULL PARAMETERS 

In order to properly determine the Weibull distribution 
which fits a given set of failure data, the values of beta and 
eta need to be estimated. This process is called Weibull 
analysis [Ref. 5] . 

Weibull analysis is best performed using Weibull graph 
paper, an example of which is in Appendix A. In performing a 
Weibull analysis, there are several questions that need to be 
answered. The primary one being, "Can the data be described 
by the Weibull distribution?" A proper Weibull plot is 
represented by a straight line, thus if the data does create 
a straight or approximately straight line (where the majority 
of the points fall in a straight line, but there may be a few 
outliers) , then the data can be represented by the Weibull 
distribution. Appendix A also contains a typical plot of 
Weibull data on Weibull graph paper. 

Weibull graph paper has unique scales for the X and Y 
axes, which graphically will create a straight line if the 
data is Weibull. The values on the X axis represent the 
measure of life for the data, and are used as the primary 
estimator for the value of eta. The Y axis represents the 
expected cumulative percentage of components that have failed. 
By measuring from the 63.2% line on the Y axis, the value of 
beta can be determined. 

Once time to failure data has been collected, it must be 
ordered from first failure to last failure. Once in the 
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propsr order, the Median rank inust be determined. Appendix B 
provides an example of a table for the determination of the 
median rank based upon the rank order and the sample size. 
Appendix B depicts the following (Sample size = 3): 

Time to Failure Median Rank 

1 = 100 hrs 20.6 

2 = 200 hrs 50.0 

3 = 300 hrs 79.3 

These values are plotted on the Weibull paper, with the 
hours on the X axis and the median rank on the Y axis. Next, 
a straight line is fitted to the data. Then eta is determined 
by drawing a vertical line down to the X axis from the point 
where the data line crosses the Y=63.2% horizontal line. This 
point on the X axis is the value of the characteristic life. 
In the case of the chart in Appendix A, that value is roughly 
240 hours. 

Beta is determined by moving down the data line to a 
point which is horizontally one inch away from the vertical 
line corresponding to the characteristic life and then 
measuring down from the intersection of the 63.2% line and the 
data line to the intersection of the horizontal and vertical 
lines. Slope is defined as rise over run, thus the distance 
in inches down to the "one-inch horizontal intersection" is 
the slope or the Beta value. In Appendix A the distance is 
1.75 inches and therefore the value of Beta is 1.75. Because 
the value of Beta is greater than one, this indicates that the 
data is from a component which in the wearout stage of the 
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bathtub curve. 


Figure 3.5 depicts the probability density 
function of this data, with a 3 = 1.75 and an r) = 240. 



Figure 3.3; Weibull PDF for Appendix A example's 
pareuneters. 

C. PROPOSED CHANGES TO THE CEAMOD PROGRAM 

The proposal to incorporate the Weibull distribution into 
the CEAMOD program has been studied by General Electric and 
reference 9 states that they feel that Weibulls are too 
complex for this model. However, a Weibull distribution has 
the capability of allowing a failure rate that will change 
with time. For a given component in a system, the Weibull 
allows either an increasing, constant, or decreasing failure 
rate. In contrast, the exponential distribution assumes only 
that the failure rate is constant. Thus, if a component has 
reached a given number of operating hours, there is no greater 
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or lesser chance of failure after this point in time than at 
any other point in its life. 

From an analytical point of view it is extremely 
difficult to incorporate Weibull generated failures into the 
CEAMOD program. The reason is that the distribution of 
failures over a given time interval is difficult to derive 
except in the case of p = 1. For that case the Poisson 
distribution can be used since times between failures are 
exponentially distributed. However, Professor James Esary of 
the Operations Research Department of the Naval Postgraduate 
School suggested that the use of simulation results should 
provide an empirical distribution which might be useful in the 
CEA model. Thus, this thesis uses a program to simulate what 
might happen if a component's times between failure were 
Weibull distributed. 

This thesis will analyze simulations of a component's 
failures over a period of the expected life of an engine. To 
provide a basis for comparison, the general operating 
parameters for this component are assumed to be those provided 
by the General Electric example (a copy of the CEAMOD output 
for that example is presented in Appendix C). One intent of 
comparing the behavior of component failures represented by 
the exponential distribution to those represented by the 
Weibull distribution, is to determine whether there is a 
significant cost difference, in the long term, based upon the 
net present value of expenditures and savings. 
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IV. SIMULATION AND SPREADSHEET MODEL DEVELOPMENT 


A. MODEL CONCEPT 

The concept for the model is to simulate failures of a 
component over its expected life. The model incorporates 
failures that occur either in an exponential manner or in a 
Weibull manner over a period of ten years. As noted in the 
previous chapter by determining the time to failure and the 
annual number of failures, it should be possible to determine 
whether the use of the Weibull distribution is feasible. 


B. SIMAN SIMULATION PROGRAM 

In order to simulate this process, a simulation program, 
written in the SIMAN language was developed. Appendix D 
contains the program used to create the failure data. 

The program begins by introducing ten engines into the 
system at time zero. It is assumed that all ten engines are 
in a "ready for flight" condition. The component of interest 
fails based upon the assumed distribution. The engine then 
enters into a repair facility, is repaired, and the engine is 
then again ready for service. That engine fails again at a 
future time based upon the assumed distribution, and the 
process is repeated. 

In order to run this program, several assumptions were 
made. These are shown in Table 4.1. The basic time 
parameters were taken from the data provided by the General 
Electric example, and others were assumed by the author. 
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PARAMETER 

VALUE 

■MfUif 

ANNUAL FLIGHT HRS 

240 

GE Example 

MEAN TIMES TO FAIL 

20,50,100,500 

GE Example 

(HOURS) 


Author 

REP FAC DELAY TIME 

0 

Author 

(HOURS) 



NBR ENGINES/YEAR 

10 

Author 

NBR OF REP FAC 

1 

Author 

LIFE OF ENG (YRS) 

10 

Author 


Table 4.1: SIMAN parameters for the failure event simulation. 

The General Electric example uses an average annual 
fight hours value of 240 hours per aircraft per year. For the 
component being analyzed in the example, the assumption is 
that the current component has a MTBF of 50 hours and the 
proposed component has a MTBF of 500 hours. In addition to 
these values, MTBF values for the exponential distribution of 
20 and 100 were chosen as well. The reason for these 
additional values was to simulate potential failures at a 
higher failure rate (20) , and those at a rate (100) between 
the two values of the example problem. The benefit of this is 
to provide additional information for failure rate 
characteristics not necessarily in line with the original 
assumptions. The eta value (characteristic life) for the 
Weibull distribution was assigned using the same values as for 
the mean of the exponential distribution. The beta factors 
were all greater than one corresponding to a component having 


24 




an increasing failure rate. The beta values ranged from 1.5 
to 5 (see Table 4.2). 

The repair facility delay time of zero hours was a 
necessary assumption to ensure the SIMAN program was 
comparable with the CEAMOD in its assumptions of repair time. 
The number of engines introduced per year, ten, and the engine 
life value, ten years, were chosen to provide sufficient data 
points to complete a credible analysis. The number of repair 
facilities, one, has no impact on this simulation, but if a 
repair delay time was introduced it could impact the number of 
failures per year by being a potential bottleneck. 

Once the parameter selection was made and the program 
written, the probability distribution of times between 
failures was assigned. As discussed above the author chose 
the following exponential and Weibull distributions to test in 
the program: 


DISTRIBUTION 

EXPO(MTBF), WEIBULL{ri, 3) 

EXPONENTIAL 

20,50,100,500 (hours) 

WEIBULL 

(20,1.5) (20,2) (20,3) (20,5) 

WEIBULL 

(50,1.5) (50,2) (50,3) (50,5) 

WEIBULL 

(100,1.5) 

(100,2) 


(100,3) 

(100,5) 

WEIBULL 

(500,1.5) 

(500,2) 


(500,3) 

(500,5) 


Table 4.2: Distributions and their parameters used for failure 
times in the simulation (all times are expressed in hours). 
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As discussed in Chapter III, the actual means and 
variances for the Weibull are different than those for the 
exponential for the same eta value. Table 5.3 reflects the 
theoretical means and variances in hours for the assumed 
distributions presented in table 4.2. These values for the 
means and variances were computed using the equations (3.4) 
and (3.5) from Chapter III. The theoretical means, for any 
value of eta, tend to decrease as beta increases until it 
reaches approximately 2.2. At this point the mean begins to 
increase slowly, asymptotically approaching the mean of the 
exponential. The variances, for all values of eta, decrease 
as beta increases. Theoretically, for a beta value of 
infinity, the mean for the Weibull equals the mean for the 
exponential and the variance is zero. 

The SIMAN program begins by creating ten "entities" 
(engines) and assigns consecutive engine numbers, "ENGNUMs", 
to each, beginning at one. The time of arrival into the 
program is assigned as the current time in the simulation. 
With the ten engines, all operationally capable, the program 
"delays" each engine by a time equivalent generated by the 
assumed distribution. This delay simulates the operational 
phase of each engine. At the end of the delay the engine 
component fails. 

Once a failure has occurred, the engine enters into a 
repair queue to await repair. Once in repair the author 
assumes that the time in the repair facility is zero. 
Although this is not a realistic assumption for real world 
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operations, it was necessary to approximate the assumptions of 
the CEAMOD program. The CEAMOD assumes that each engine fails 
according to the exponential distribution, and does not assume 
that the repair channel may be a delaying factor in the return 
to operational capability. Thus the author has chosen a delay 
of zero for the repair channel. This SIMAN program can, of 
course, be easily modified to reflect the reality of repair 
channel delays. 

After the repaired engine exits the repair channel, the 
time to failure is tallied for each engine. The engine is 
then returned to the point where the next delay time is 
generated. It is then delayed again by that amount and the 
process repeats. 

Based upon the assumptions stated in Table 4.1, the 
program runs for a total of 2400 hours for each engine; that 
is, 240 hours per year for ten years. Because there is no 
"delay" for each engine in the repair channel, the "system" is 
in a steady state with no need for an initial warm-up period. 
If a delay time for repair were introduced, the SIMAN program 
would need to include some time for a warm-up in order to 
achieve steady state. 

The data that was written to the tally files reflects the 
point in time when each failure occurred. An example output 
from the program is included in Appendix E. This output 
provides the program user with several important pieces of 
data. For each engine it lists the average time between 
failures (MTBF), the variation of the simulated times, the 
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Distribution 


MTBF (Hours) 


Variance (Hours) 


Expo(20) 


20.0 


400.0 


Weib(20,1.5! 


18.1 


150.3 


Weib(20,2) 


17.7 


85.8 


Weib(20,3! 


17.9 


42.1 


Weib(20,5) 


18.4 


17.7 


Expo(50) 


50.0 


2500.0 


Weib(50,1.5! 


Weib(50,2! 


45.1 


44.3 


939.2 


536.5 


Weib(50,3) 


44.6 


263.3 


Weib (50,5) 


45.9 


110.6 


Expo(100! 


Weib(100,1.5) 


100.0 


90.3 


10000.0 


3756.9 


Weib(100,2! 


88.6 


2146.0 


Weib(100,3) 


89.3 


1053.3 


WeibdOO, 5) 


Expo(500) 


91.4 


500.0 


442.3 


250000.0 


Weib(500,1.5) 


451.4 


93922.6 


Weib(500,2) 


443.1 


53650.5 


Weib(500,3) 


Weib(500,5) 


446.5 


459.1 


26333.2 


11057.5 


Table 4.3; Theoretical means and variances for exponential 
and Weibull distributions. 
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minimum, maximum, and the total numher of failures that could 
be expected over the life of 2400 hours. 

The progratn also writes the actual failure times to dat 
files, as specified in the experimental frame. These "dat" 
files are important in analyzing the behavior of the failures, 
and the effect of the assumed statistical distribution. An 
example of the output from one file is included in Appendix F. 
The "dat" files can be converted to a spreadsheet readable 
format by using the "export" feature of SIMAN. 

After importing the data into a spreadsheet, it can be 
easily analyzed and graphed to reflect the potential number of 
failures each year over the life of the engines. By reading 
the data from all ten engines into the spreadsheet, and 
grouping the number of failures by years (group together the 
failures from all engines, based upon the time of failure) it 
is easy to determine how many engines are expected to fail in 
any given year. These groupings can then be graphed to 
further visualize the number of failures each year over the 
life-cycle. Examples of such graphs are included in Appendix 
G. 

The use of a simulation program is one way to generate a 
distribution of failures, and not necessarily the easiest way. 
By using a spreadsheet, such as Excel, a distribution can be 
created for the Weibull distribution. The benefit of this 
method is that it can be easily accomplished without the 
complication or expense of a simulation program. The drawback 
is that other parameters, such as a repair channel delay time 
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cannot be introduced. 

The process for creating the Weibull failure times can be 
accomplished with the following equation:. 

Weibulhx] . 

The value (-In(rand)) is a spreadsheet command that takes the 
negative of the natural log of a randomly generated, uniformly 
distributed number between zero and one. By copying this 
equation down a column in the spreadsheet, an entire set of 
failures can be created. This equation creates a time to 
failure for each cell. Column D of table 4.4 is an example of 
the output from this process. In that table column E is then 
used to show the cumulative clock time. This is necessary to 
determine the cut-off point for failures over a specific 
period of time. Counting up the number of failures from 
Column D provides the number of failures over a given period 
of time. From Table 4.4 it is easy to see that during a 
period of 240 hours, there were seven failures. That is 
determined by simply counting the number of failures that 
occurred prior to the "cut-off" time of 240. 

In order to analyze many engines over many years, the 
equation for generating failures can be copied to multiple 
columns, one for each engine. To determine the distribution 
of number of failures over all engines per year, the failure 
time columns can be analyzed using the "data analysis" 
function in Excel, and histograms representing number of 
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failures in a year can be plotted. This process is quite 
easy, and by changing the values of beta and eta, the 
distributions for all Weibull parameters can be determined. 


A 

B 

C 

D 

E 

F 

1 

BETA = 1 


Failure 

Cum Time 


2 

ETA = 50 


28.97 

28.97 


3 



19.77 

48.74 


4 



14.33 

63.07 


5 



9.76 

72.83 


6 



54.04 

126.87 


7 



32.16 

159.03 


8 



62.51, 

221.54 


9 



58.78 

280.32 



Table 4.4; Example of spreadsheet generated failure times. 


In order to determine the potential effects of a given 
distribution, relative to the current CEAMOD program's 
assumptions, the data must be analyzed from a financial 
perspective. That is where a modified, simplified spreadsheet 
analysis comes into play. The next section will address this 
function. 

C. EXCEL FINANCIAL SPREADSHEET PROGRAM 

The financial spreadsheet developed by the author 
computes the net present value (NPV) of the proposed changes 
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over the life of the program. Appendix I contains examples 
of the spreadsheet program. The following are the specific 
assumptions and calculations included in the program. 

1. ASSUMPTIONS AND CRITICAL VALXJES (TOP LEFT CORNER OF 
SPREADSHEET) 

a. Interest Rate (Row 1) 

The assumed interest rate is 10% and is the discount 
rate at which the cost or savings from the modification will 
be valued for all future years. Thus a cost of $1000.00 next 
year is equivalent to $909.09 today. This is obtained by 
dividing the cost of $1000.00 by (1 + interest rate). For all 
future years, (1 + interest rate) is raised to the power of 
the number of years. Thus the value of $1000.00 in five years 
is equivalent to 1000/(1.10)® or $620.92 today. 

b. Flight Hours per Year (Row 2) 

This value is defined as the expected number of 
flight hours each engine will operate for in each calendar 
year. Thus, for example, with ten engines operating, that 
value is multiplied by ten. 

c. Cost of Modification (Row 3) 

This cost is the expected cost of the modified 
component which is in excess of the current component's 
replacement cost. This value can be positive or negative 
depending upon the specific component. This value does not 
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includs r©s 6 a. 3 rch and devslopitisnt costs. It is simply th.6 
(jiff 0 ]^ 0 ];iQ 0 p 0 twssn ths cost of a nsw unit of ths old componsnt 
and the cost of a new unit of the new component. 

d. Labor Cost (Row 4) 

The cost per hour to have trained technicians 
replace the failed component. This price is an aggregate for 

all technicians needed to perform the job. This labor cost 

/ 

may be at either the organizational level, intermediate level, 
or depot level or all three. The assumption of $25.00 per 
hour is the organizational level labor rate. 

e. Current Failure Rate (Row 5) 

This rate is the expected number of failures per 
engine per flight hour. It is equivalent to one divided by 
the Mean Time Between Failures in flight hours for the 
exponential distribution. Therefore if the failure rate is 
.02 per flight hour then the MTBF is 50 flight hours. With 
240 flight hours per year, the expected number of failures is 
(.02 * 240) or 4.8 failures per year assuming the exponential 
distribution for times between failures. The value displayed 
in Appendix I is the expected number of failures per flight 
hour, rounded to three significant digits. 


f. Proposed Failure Rate (Row 6) 

The proposed failure rate is defined in the same manner 
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as the current failure rate. It represents the failure rate 
of the proposed design change. This value is used to 
determine the total expected failures of the proposed 
configuration per year, and is rounded to three significant 

digits. 


g. LaJaor Hours (Row 7) 

This represents the number of hours of labor which will 
reguired to repair an engine once it has failed. The labor 
hours multiplied by the Labor cost per hour represents the 
total cost (of labor) to repair each engine. For the example 
in Appendix I, ten is the assumed number of labor hours. This 
value could be incorporated into the SIMAN model as the repair 
channel delay time. 

h. Current Failures per Year (Row 8) 

This value represents the total number of expected 

failures per year per engine based upon the current 

configuration. This value is obtained from the output of the 
SIMAN program by taking the average number of failures for all 
engines over the ten-year period. If there are 1000 failures 
for ten engines in ten years, that equates to an average of 
ten failures per engine per year. 

i. Projected Failures per Year (Row 9) 

This value represents the total number of expected 

failures per year per engine based upon the proposed 
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QQnfiguifation. This vslus is obtsinsd. firotn ths output of tho 
SIMAN program and is determined in the same manner as the 
current failures per year. 

2. Spreadsheet Calculations 

a. Year (Column 1) 

This is the designation of the year in the life-cycle 
analysis: it begins at year one and continues through year 
twenty. Based upon the assumed engine life, an engine 
introduced in year one will be phased out of service ten years 
later, in year eleven. 

b. Engines Introduced (Column 2) 

This represents the number of engines with the improved 
component introduced in a given year. For sake of simplicity 
the author has chosen to introduce ten engines in each of the 
first ten years of the cycle. Each engine component is 
assumed to be introduced at the beginning of the year. 

c. Attrition (Column 3) 

This represents the number of engines or engines with 
specific components removed from service in a given year. The 
author has assumed a ten-year engine life and therefore the 
calculations show, for example, the engines being introduced 
in year one are being removed at the start of year eleven. 
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d. Total in Service (Column 4) 

This represents the total number of engines operating 
during a given year. It is the sutn of all previous 
introductions minus all previous attritions. All 
introductions and attritions are assumed to occur at the start 
of the year. 


e. Total Flying Hours (Column 5) 

This represents the expected total number of flying hours 
performed in a given year for all operational aircraft. This 
number is obtained by multiplying the expected number of 
flying hours per aircraft per year by the total number of 
engines in service during that year. 

f. Current Expected N\nnber of Failures (Column 6) 

This represents the total number of failures that can be 
expected in a given year for the current configuration. This 
value is obtained by multiplying the total flying hours value 
by the current failure rate (row 8 of the assumption section). 
It may appear that the value in column 6 does not match this 
calculation because the value displayed for the current 
failure rate is only displayed to three significant digits. 
The actual value used is the complete value of the current 
failure rate (ten digits) without rounding. 

g. Proposed Expected Number of Failures (Column 7) 

This value represents the total number of failures that 
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can be expected in a given year for the proposed 
^on. This value is obtained in the satue manner as 
column 6 except that the "proposed" value (row 9 of the 
assumptions section) is used vice the "current" value (row 8 
of the assumptions section). 

h. Current Costs (Column 8) 

This value represents the current year dollar value of 
expected outlays for maintenance of the current configuration 
due to failures. This figure is obtained by multiplying the 
Current Expected Failure values (column 6) by the labor cost 
and the labor hours value. In the first example in Appendix 
I, for year one, the current cost equals (120.7 failures x $25 
per hour x 10 hours) or $30,175.00. The sum at the bottom of 
the column represents the total costs over the entire life 
cycle of operating an engine without a component modification. 

i. Projected Costs (Column 9) 

This value represents the current year dollar value of 
the expected outlays for maintenance due to failures of the 
proposed configuration. It is obtained by multiplying the 
number of engines introduced (column 2) by the cost of a 
modification, and adding that to the projected expected 
failures (column 7) multiplied by the product of labor cost 
and labor hours. The entry in the first example in Appendix 
I for year one equals 

(10 new engines x $1000 mod cost) + (47.2 expected failures x 
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$25.00 labor cost x 10 labor hours) which totals $21,800.00. 
The sum at the bottom of the column represents the total 
projected costs over the entire life cycle of operating an 
0 ngine with a component modification. 


j. Projected Savings (Column 10) 

This value is the difference between the Current Costs 
and the Projected Costs (column 8 minus column 9). This value 
may be either positive or negative depending upon the overall 
cost differential. The sum at the bottom of this column 
represents the total savings over the entire life cycle of 
operating engines with the modification instead of without the 
modification. 


k. NPV of Projected Savings (Column 11) 

This value is the cumulative discounted value of the 
projected savings in each year. Since all costs are assumed 
to occur at the beginning of the year, the first year's value 
is not discounted. The present value of second year's cost 
represents the value of the projected savings for year two 
divided by (1 + interest rate)y^^^-\ plus the savings from the 
previous year. In the first example in Appendix I, the 
savings for year three equals (45,125/(1.1)^ +$32,693.00 = 
$69,986.57. The value at the bottom of the column represents 
the total net present value of the savings from making the 
configuration change. This is usually the key value of the 
CEA model. This value is used by the decision maker in 
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deciding whether it is worthwhile to make the change. 


D. SUMMARY 

By using the SIMAN program developed by the author, and 
the EXCEL spreadsheet program, it is possible to determine 
whether there can be any expected differences between the use 
of the Exponential and Weibull distributions. 

The simulation of the life cycle of engine failures is 
critical to this analysis. By looking at the simulated 
failures of engines based upon an assumed distribution, it is 
possible to project the future. In an operational arena, 
through the use of Weibull analysis, the engineer would, with 
reasonable certainty, know what the failure distribution would 
be. By then using this distribution in a simulation program, 
he/she could see what would be expected in the long run. 

The CEAMOD program has assumed that the distribution is 
exponential, even though the Weibull may be a more logical 
choice. Components of mechanical equipment do tend to wear 
out as time goes by. The Weibull allows the analyst to 
reasonably determine what is the expected life remaining of a 
component once it has operated for a given number of hours. 

In the next chapter, the output data from the simulation 
program and the financial spreadsheet will be analyzed, and a 
comparison will be made considering the various parameter 
changes. For example, by comparing the expected number of 
failures from the exponential and the Weibull distributions, 
over the life cycle of the engines, it can be determined 
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whether there is any merit to changing to a Weibull-based 
analysis. 
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V. ANALYSIS OF OUTPUT 


A. INTRODUCTION 

This chapter will focus on the output provided by the 
SIMAN and financial spreadsheet programs. The purpose of this 
section is to compare the output data from the simulation 
using the CEAMOD assumption of an exponential distribution 
versus a Weibull distribution. 

B. SIMAN PROGRAM OUTPUT 

In order to produce the necessary output data for each of 
the relevant distributions (see Table 4.2), it was necessary 
to update the program after each run. The only change which 
needed to be made was to modify the "Delay" factor to reflect 
the distribution needed to produce the results. 

The first analysis was to look at the simulated 
MTBFs of the components over the expected life of the engines. 
In order to determine the validity of the simulation model it 
was necessary to compare the SIMAN results to the theoretical 
means and variances. Table 5.1 provides that comparison. 
Based upon the comparison of the means and variances from a 
simulation of 100 engine years (ten engines for ten years 
each) to the means and variances of the theoretical 
distributions, it appears that the simulation replicated 
theory very closely. The highest difference between means 
occurred for Weib(500,2) but the simulated mean was only 5.7% 
lower than the theoretical mean. The highest variance 
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E (20) 


W(20,1.5) 


W(20,3) 


W(20,5) 


E (50) 


W(50,1.5) 


W(50,2) 


W(50,3) 


W(50,5) 


E (100) 


W(100,1.5) 


W(100,2) 


W(100,3) 


W(100,5) 


E (500) 


W(500,2) 
W(500,3) 
W(500,5) 


OBS 

MEAN 


19.8 


18.0 


17.9 


18.4 


47.1 


45.5 


44.3 


44.9 


46.1 


108.0 


95.3 


92.0 


91.6 


93.3 


517.2 


419.1 

447.8 

458.8 


THEOR 

DIFF 

OBS 

THEOR 

MEAN 

(0-T) 

VAR 

VAR 

20.0 

( .2) 

355.7 

400 

18.1 

.1 

143.0 

150.3 

17.7 

0 

82.7 

85.8 

17.9 

0 

40.4 

42.1 

18.4 

0 

16.6 

17.7 

50.0 

(2.9) 

1964.4 

2500 

45.1 

.4 

890.9 

939.2 

44.3 

0 

521.6 

536.5 

44.6 

.3 

259.8 

263.3 

45.9 

.2 

108.3 

110.6 

100.0 

8 

9967.3 

10000 

90.3 

5 

3658.3 

3756.9 

88.6 

3.4 

2056.3 

2146.0 

89.3 

2.3 

1028.7 

1053.3 

91.4 

1.9 

419.4 

442.3 

500.0 

17.2 

155697 

250000 

451.4 

(14.4) 

71536.1 

93922.6 

443.1 

(24) 

43727.4 

53650.5 

446.5 

1.3 

22772.4 

26333.2 

459.1 

( .3) 

8588.0 

11057.5 


DIFF 

(0-T) 


(44.3) 


(7.3) 


(3.1) 


(1.7) 


( 1 . 1 ) 


(535.6) 


(14.9) 


3.5 


(2.3) 


(32.7) 


98 


(89.7) 


(24.6) 


(22.9) 


(94303) 



(9923.1) 
(3560.8) 
(2469.5) 


Table 5.1: Comparison of observed to theoretical means and 
variances. 
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difference occurred for Expo(500) and the simulated variance 
was 60% lower than the theoretical variance. Both differences 
are mainly due to the fact that the mean times between 
failures are quite large and therefore longer simulation runs 
will be necessary before we could expect the simulated means 
and variances to converge to the theoretical values. 

The second analysis undertaken was to determine whether 
the number of failures cycled over the years of the engine 
life cycle. The SIMAN program was run multiple times to 
simulate the potential number of times engines failed over 
each year. Appendix G provides a graphical illustration of 
the output from this process. Table 5.2 adds some additional 
information about the results. 

As is graphically depicted in Appendix G, the "curve" of 
failures has the greatest "cyclic" nature when the exponential 
distribution is assumed. When the Weibull is used, even 
assuming the same characteristic life parameter, and the value 
of beta increases, the number of simulated failures in each 
year tends to smooth out. 

This can best be interpreted by saying that for 
components that have been determined to be "wearout" 
candidates, and after a Weibull analysis has been conducted to 
determine the appropriate parameters for the Weibull 
distribution, the higher the beta value, the more likely it is 
that there will be a more constant number of failures per 
year. In the case of an assumption of a high characteristic 
life, 500 hours, the smoothing effect is less pronounced, but 
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DISTRIBUTION 

AVG NBR OF 

COMPONENT 

FAILURES/YR 

— 

MAX NUMBER OF 

FAILURES IN 

ANY YEAR 

MIN NUMBER OF 

FAILURES IN 

ANY YEAR 

EXPO(20) 

120.7 

144 

102 

WEIB(20,1.5) 

133.1 

144 

124 

WEIB(20,2) 

135.3 

144 

128 

WEIB (20,3) 

133.9 

139 

126 

WEIB(20,5) 

129.8 

133 

124 

EXPO (50) 

47.2 

62 

35 

WEIB(50,1.5) 

52.4 

59 

48 

WEIB (50,2) 

53.4 

58 

50 

WEIB(50,3) 

53.1 

59 

50 

WEIB(50,5) 

51.7 

54 

48 

EXPO(100) 

22.0 

27 

11 

WEIBdOO, 1.5) 

25.1 

30 

20 

WEIB(100,2) 

25.8 

29 

22 

WEIBdOO, 3) 

25.1 

29 

22 

WEIBdOO, 5) 

25.3 

27 

21 

EXPO(500) 

5.3 

10 

0 

WEIB(500,1.5) 

5.2 

8 

2 

WEIB(500,2) 

5.3 

8 

3 

WEIB(500,3) 

5.1 

7 

1 

WEIB(500,5) 

4.9 

8 

0 


Table 5.2: S-ummary of simulated engine failures. 
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is still present. For components with a low MTBF (a high 
failure rate) the presence of the smoothing trend is quite 
pronounced. 

From the data generated by the SIMAN simulation program, 
ttie author concludes that the Weibull distribution does not 
show an increase in the "cyclic" nature of the numbers of 
failure data, but, in fact shows a decrease in that 
phenomenon. The Weibull distribution for beta values greater 
than one appears to have a "smoothing" effect on the number of 
failures. 

The third part of the analysis was to determine the 
nature of the distribution of the annual number of failures 
when the time between failures is Weibull distributed. 
Appendix H presents histograms which can be used to identify 
some characteristic distributions. The data for these plots 
was developed using the simulation program, and they represent 
the number of failures occurring each year for 100 engine 
years. The means and variances for each set of plots are 
summarized in Table 5.3. This table shows the mean number of 
failures per year for the Weibull decreases initially as the 
beta value increases. However it does begin to increase again 
once the beta value exceeds approximately 2.2. This 
phenomenon occurs for the same reason that the MTBF decreased 
until beta reached approximately 2.2 and then it increased as 
shown in Table 4.3. The variance of the number of failures, 
however, is a strongly decreasing function of beta. 

This table also shows that the means and variances for 
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DISTRIBUTION 

MEAN NBR OF 

FAILURES/YEAR 

VARIANCE OF 

FAILURES/YEAR 

EXPO(20) 

12.07 

10.22 

WEIB(20,1.5) 

13.31 

5.02 

WEIB(20,2) 

13.53 

2.87 

WEIB(20,3) 

13.39 

1.83 

WEIB(20,5) 

12.98 

.80 

EXPO(50) 

4.72 

4.14 

WEIB(50,1.5) 

5.24 

2.34 

WEIB(50,2) 

5.34 

- 1.09 

WEIB(50,3) 

5.32 

.94 

WEIB(50,5) 

5.17 

.46 

EXPO(IOO) 

2.2 

1.95 

WEIB(100,1.5) 

2.51 

1.32 

WEIB (100,2) 

2.58 

.83 

WEIB(100,3) 

2.51 

.47 

WEIB(100,5) 

2.53 

.29 

EXPO (500) 

.53 

.49 

WEIB(500,1.5) 

.52 

.35 

WEIB(500,2) 

.53 

.31 

WEIB(500,3) 

.51 

.25 

WEIB(500,5) 

.49 

.27 


Table 5.3: Mean number of failures per year and variance 
of nxamber of failures per year. 
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the exponentially distributed failures are close in value but 
not equal. They are theoretically equal since the number of 
failures per year is Poisson. The slight differences here are 
consequences of a limited number of simulation runs. 

On the first page in Appendix H, for those graphs with an 
eta value of 20, the exponential distribution represents an 
approximate Poisson distribution. In the case of each of the 
"Weibull plots" they do not look like a Poisson distribution 
or even a discrete form of a Normal distribution when beta is 
less than two. For beta values greater than two a Normal 
distribution using a continuity correction might be a 
reasonable approximation. 

To better discern the nature of the transition from the 
Poisson for an eta of 20 and the approximate Normal for beta 
greater than two, two additional plots have been included on 
the second page of Appendix H. These plots have beta values 
of 1.25 and 1.75, and an eta value of 20. The mean number of 
failures per year for the Weibull (20,1.25) is 12.87 with a 
variance of 7.12. For the Weibull(20,1.75) the mean number of 
failures is 13.48 with a variance of 5.24. 

The histograms for eta of 50, 100, and 500 are presented 
on the remainder of the pages in Appendix H and show an 
increasingly better Poisson distribution for the exponential 
failure times than was the case for the eta of 20. The 
transition to the Weibull and the influence of beta is also 
increasingly "smoother". 

The fourth aspect of this analysis is to determine 
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wh6th6ir ths Wsibull Tna.y hsvs sny impsct on financis.! docision 
making. That will be addressed in the next section. 

C. FINANCIAL SPREADSHEET ANALYSIS 

One of the principal goals of the CIP program is to 
maintain an engine design which allows the maximum aircraft 
availability at the lowest cost to the government [Ref. 10] . 
Thus the question remains, "does the introduction of the 
Weibull distribution into the CEAMOD program more accurately 
the reflect the costs and/or potential savings?" By 
introducing the expected failure rates per year into the 
financial spreadsheet developed by the author, that question 
can be answered. 

In order to most closely reflect the General Electric 
example, shown in Appendix C, the author has chosen to use the 
exponential MTBF of 50 as the base, and compare the other 
three values (20, 100, and 500) to that base. The only 
exception is the assumption that in making a comparison the 
designers would desire to move from a lower MTBF to a higher 
MTBF. Thus, in the comparison of 20 hours, that will be the 
"current" and 50 will be the "proposed". In the other two 
comparisons 50 is the "current" and either 100 or 500 is the 
"proposed". 

The author has also chosen to first compare similar 
distributions. Thus the exponentials will be compared and 
the Weibulls with like beta values will only be compared. 
After that a comparison will be made between the exponential 
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and Weibull for unlike parameters. In that comparison 
exponential will be chosen as the "current" distribution and 
one of several Weibull distributions will be chosen for the 
"proposed". 

The highlighted rows in Table 5.4 indicate the highest 
projected saving for each set of distribution comparisons. 
Example spreadsheets from which this data was extracted are 
provided in Appendix I. The savings for a Weibull 
distribution follow a similar pattern to the number of 
failures per year, increasing as beta increases to two, then 
slowly decreasing until they reach the level achieved by the 
exponential. 

Based upon the projected costs/savings table, it is clear 
that given all assuitptions previously discussed, there can be 
a significant savings differential between the projections 
based on both the exponential distribution and the Weibull. 
However, the savings between Weibulls is larger because the 
difference between the expected number of failures per year 
for the current configuration and the proposed configuration 
is larger with the Weibull assxomption. Thus, one can expect 
a greater decrease in the number of failures per year with a 
Weibull assuit^tion for both the current and proposed 
configurations. This greater decrease translates into an 
increase in potential savings. For example, the difference 
between the expected annual number of failures for Expo(20) 
and Expo (50) is 7.35 (12.07 for Expo (20) minus 4.72 for 
Expo(50)). For the comparison between Weib(20,1.5) and 
Weib(50,1.5) the difference is 8.07 (13.31 for Weib(20,1.5) 
minus 5.24 for Weib(50,1.5)). 
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CURRENT 

PROPOSED 

PROJECTED 

FIRST YEAR OF 

DISTRIBUTION 

DISTRIBUTION 

SAVINGS 

POS CASH FLOW 

EXPO (20) 

EXPO (50) 


1 

WEIB (20,1.5) 

WEIB (50,1.5) 

$854,092 

1 


WEIB (20,2) 


WEIB (20,3) 


WEIB (20,5) 


EXPO (50) 


WEIB (50,1.5) 


WEIB(50,2) 


V«EIB(50,3) 


WEIB(50,5) 


EXPO(50) 


WEIB(50,1.5) 


WEIB(50,2) 


V1EIB(50,3) 


WEIB(50,5) 


WEIB (50,2) 


WEIB (50,3) 


WEIB (50,5) 


EXPO (100) 


WEIB(100,1.5) 


WEIB(100,2) 


WEIB(100,3) 


WEIB(100,5) 


EXPO(500) 


WEIB(500,1.5) 


WEIB(500,2) 


WEIB(500,3) 


WEIB(500,5) 


$867,797 


$854,092 


$824,397 


$ 220,221 


$244,205 


$247,632 


$253,342 


$233,926 


$410,953 


$471,485 


$481,764 


$481,764 


$466,917 



Table 5.4: Projected savings from changes in distribution 
parameters. 

Table 5.5 provides a comparison of the projected 
costs/savings when changing from the exponential distribution 
to the Weibull. This might correspond to a situation where the 
old component had a constant failure rate and the new one did 
not. The values for the first two lines are negative because 
the comparison is between distributions with the same eta 
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CURRENT 

DISTRIBUTION 

PROPOSED 

DISTRIBUTION 

PROJECTED 

(COSTS)/SAVINGS 

EXPO (20) 

WEIB(20,1.5) 

(209,211) 

EXPO (20) 

WEIB(20,5) 

(171,522) 

EXPO (20) 

WEIB(50,1.5) 

712,470 

EXPO (20) 

WEIB(50,2) 

701,049 

EXPO (50) 

WEIBdOO, 1.5) 

184,816 


Table 5.5; Expected cost/savings conparison from changing 
from an exponential to a Weibull. 

value. In these cases the average number of failures per year 
increases, thus there is an additional cost for this change. 
However, the cost is less for a WEIB(20,5) than for a 
WEIB(20,1.5) because the niomber of expected failures per year 
is less for the WEIB(20,5). For the next two lines, the 
comparison between an EXPO(20) and a WEIB(50,1.5) and 
WEIB(50,2) shows a savings, but not as great a savings as the 
change shown in Table 5.4 when the change was made with like 
distributions. This is due to the increase in the expected 
number of failures per year between distributions. The 
difference between the average number of failures for a 
comparison of EXPO(20) and EXPO(50) is 7.35. However, the 
difference between EXPO (20) and Weibull(50,1.5) is 6.83 and 
for Weibull(50,2) it is 6.73. This decrease in the difference 
between the current and proposed distributions accounts for 
the decrease in the projected savings. The savings are 
primarily achieved from the shift from an MTBF of 20 hours to 
one approaching 50 hours vice a shift to a different 
distribution. The same is true for the last comparison in 
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Table 5.5. 

Based upon this analysis, the author concludes that there 
is a significant difference in costs between the exponential 
and Weibull distributions. In addition, when applied to 
Weibull with like beta values but differing eta values, those 
savings can be quite large. The savings for a Weibull 
distribution follow a similar pattern to the number of 
failures per year, increasing as beta increases to two, then 
slowly decreasing until they reach the level achieved by the 
exponential. 
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VI. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

The key questions that needed to be answered by this 

thesis research effort were: 

1. Does the assumption of the Weibull distribution mean 
that there will be "cyclic" numbers of failures over the years 
of the engine life-cycle? 

2. What does the probability distribution for the number 
of failures in a year look like if failure times are Weibull? 

3. Is there any significant difference in the net 
present value of future life cycle costs/savings by assuming 
the Weibull instead of the exponential distribution for the 
distribution of times between failures for an aircraft engine 
component? 

The exponential is a special case of the Weibull and it 
is well known that if times between failures are exponentially 
distributed then the number of failures over a given period of 
time are Poisson distributed. However, for any other case of 
the Weibull the distribution of the number of failures over a 
given period of time cannot be derived analytically. 
Therefore, the use of simulation is needed to study the nature 
of this distribution. As a consequence, the author was able 
to generate the empirical distributions in question for a 
range of Weibull parameters. And, in addition, obtain an 
estimate of the mean number of component failures per year for 
different parameter values of the Weibull. That information 
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was then used in a simple version of the CEA model to compare 
the results derived from assuming exponentially distributed 
times between failures with the results from assuming Weibull 
distributed times between failures. 


B. CONCLUSIONS 

From the results of this analysis the three questions 
listed above can be answered. The answer to the first 
question regarding "cyclic" numbers of failures is that the 
use of the Weibull, with a shape factor greater than one, 
leads to less variation in the annual simulated number of 
failures than the exponential. As the beta value (shape 
factor) increases the variance in the total number of failures 
is reduced. However, the Weibull assumption with a shape 
factor greater than one does increase the expected number of 
failures per year since the mean time between failures becomes 
smaller. 

With the data analyzed for this thesis, the distribution 
of the number of failures per year for the Weibull 
distribution appears to approximate the Normal distribution 
with beta values greater than two. When the beta value is one 
(the exponential distribution) the Poisson distribution is 
clearly represented and can be theoretically justified. With 
beta values between one and two, the distribution appears to 
be transitioning from Poisson to a Normal approximation. As 
the beta value increases the skewness of the distribution, and 
the variance decrease and it begins to approach an approximate 
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Normal distribution although this remains to be statistically 
tested. 

The question regarding differences in the net present 
value of the life cycle costs of an aircraft engine under a 
proposed component change also produces a positive result. 
If, in fact, the component is failing in a manner that 
reflects a Weibull distribution and, there is a proposal to 
change a component to one that has a higher MTBF, the use of 
the Weibull may show a greater net present value of the 
savings than if the change was from assuming the component 
failed according to the exponential distribution. Based upon 
the analysis in this thesis those savings may be significant. 


C. RECOMMENDATIONS 

The author recommends that the use of the Weibull in the 
CEAMOD program continue to be evaluated. The first step in 
this process is to determine if a theoretical distribution can 
describe the number of failures per year. A Chi-square 
goodness of fit test of potential distributions should be 
made. Next, it is appropriate to determine what kinds of 
components actually do display the characteristics of the 
Weibull. One way is through Weibull analysis as shown in 
Appendix A. If a component has its time between failures 
described by a Weibull then incorporation into the model is 
justified, even if the number of failures distribution cannot 
be determined precisely. 

In order to incorporate the Weibull into the CEAMOD, the 
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progratu designers will need to rewrite portions of the input 
and calculation sections of the model. A failure rate can not 
be used as an input value as is currently done for the 
exponential. The expected number of failures in a given time 
period will need to be the input value if the failures are 
Weibull. This can be estimated either by additional 
simulation runs, or by using a spreadsheet generated Weibull 
distribution with known beta and eta values (determined 
through Weibull analysis) as described in Chapter IV. Once 
this new input value has been determined, the calculations in 
the financial portion of the model can be modified to 
accommodate the new input data. 

Finally, this thesis has also considered the variance of 
the exponential and Weibull distributions. It has shown that 
increased beta values (greater than one) result in reduced 
variance of the number of failures over a given period of 
time. This suggests that an analysis of the distribution of 
life cycle costs may be the next step in the evolution of the 
CEA model. 
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appendix A; WEIBULL PAPER AND GRAPH 


This first page of this appendix contains an example of 
blank Weibull graph paper taken from Reference 5. This paper is 
used to plot actual failure data to determine whether the data 
characteristics fit the Weibull distribution. 

The second page contains a sample plot of three points, 
which is full explained in Chapter III. 
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The median rank table assigns the median value to each item 
in a sample size. The median rank is used for plotting against 
th X axis of Weibull graph paper. The highlighted value 3, is used 
as an example in Chapter 111. 


10 

6.6 

16.2 

25.8 
35.5 

45.1 

54.8 
64 4 

74.1 
83.7 
93.3 


20 

3.4 
8.2 

13.1 
18.0 
22.9 

27.8 

32.7 

37.7 

42.6 
47.5 

52.4 
57.3 

62.2 
67.2 
72.1 
77.0 

81.9 

86.8 

91.7 

96.5 
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APPENDIX C: SAMPLE CEAMOD OUTPUT REPORT 


This appendix contains the standard three page output report 
from the CEAMOD program. The input data is taken from Reference 
4, and is the example data provided by General Electric. 
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TITLE 

ENGINE MODEL Fwwv-xx-yyy 
TASK/ECP 


F-zz 


CEA VERSION 2 1 12/9/94 

Pg 1 


Tssk /ncorpor^tfon Input 

1 0 Incorporation Style (1.2 or 3) 

1 = Attrition 

2 = Retrofit at 1st Opportunity 

3 = Forced Retrofit Kits / Month 


2 0 Does Kit Cost Replace Normal Matnt Material Cost? 1=Yes 0=No 


3 0 Delta Production Cost 

4 0 Kit Hardware Cost - $ / Engine 

5 0 Kit Labor Manhours at O&l 

6 0 Kit Labor Manhours at Depot 

7 0 Technical Pubs Cost - Total $ 

8 0 TCTO Cost - Total $ 

9 0 Tooling/Support Equipment Cost-Total $ 

10 0 Spare Parts Factor 

11 0 Scheduled % Events being Modified 

12 0 Unscheduled % Events being Modified 

13 0 Unscheduled Event Rate allowing Modification 

14 0 Production Incorporation Date 

15 0 Field Incorporation Date 


$10,000 

$15,000 

2 

20 

$500 

$1,500 

$500 

0 % 

100 % 

100 % 

0020 

Year -^ 1993 

Year -^ 1993 


Scheduled Input 

16 0 Scheduled Maintenance Interval (TAC's) 

17 0 Calculated Scheduled Maintenance Interval Rate/1000 EFH 

18 0 Scheduled Manhours at O level 

19 0 Scheduled % Removed at O&l level 

20 0 Scheduled Manhours to Remove/Repiace at O level 

21 0 Scheduled Manhours at I level 

22 0 Scheduled % at O&l requiring Repair 

23 0 Scheduled Repair Cost at O&l level 

24 0 Scheduled % Returned to Depot 

25 0 Scheduled Manhours at Depot 

26 0 Scheduled % at Depot requiring Repair 

27 0 Scheduled Repair Cost at Depot 

28 0 Scheduled % Scrapped 

29 0 Hardware Cost to Scrap 

30 0 Scheduled Engine Test Time 
Unscheduled Input 

310 Unscheduled Event Rate/1000 EFH 

32 0 Unscheduled Manhours at O level 

33 0 Unscheduled % Removed at O&l level 

34 0 Unscheduled Manhours to Remove/Replace at O level 

35 0 Unscheduled Manhours at I level 

36 0 Unscheduled % at O&l requiring Repair 

37 0 Unscheduled Repair cost at O&l level 

38 0 Unscheduled % Returned to Depot 

39 0 Unscheduled Manhours at Depot 

40 0 Unscheduled % at Depot requiring Repair 
410 Unscheduled Repair Cost at Depot 

42 0 Unscheduled % Scrapped 

43 0 Hardware Cost to Scrap 

44 0 Unscheduled Engine Test Time 

45 0 Unscheduled Secondary Damage Costs 

46 0 Unscheduled Incidental Costs 

47 0 Number of P/N’s 
Optional Input 

48 0 % Improvement in Specific Fuel Consumption from Current to Proposed 

49 0 Aircraft Loss Rate Improvement / 1,000,000 EFH 


f Standard Inputs ] 

Fiscal Year Dollars 


1991 

NPV Rate 


10% 

Labor Cost / Manhour at O&l 

$32 32 

Labor Cost / Manhour at Depot 

$43.30 

Cost to introduce new P/N - $ / PN 

$1,524 

Cost to Maintain each P/N / Year 

$250 

Fuel Cost / Gallon 


$0 61 

Test Fuel - Gallons / Hour 

150 

Flight Fuel - Gallons / Hour 

150 

EFH / Year 


240 

TAC / EFH Ratio 


300 

TOT / EFH Ratio 


1 50 

Aircraft Cost 


$0 

Month _^ 5 


Month _^ 8 


; CURRENT I 

PROPOSED 


I 3000 ; 

4000 

T 

I 1.000 I 

0 750 


1 00 I 

00 


] 100% I 

100% 


I 10.0 ' 

100 


1 25,0 1 

25.0 


1 100% I 

100% 


I $500 I 

$500 


I 100% 1 

100% 


1 10,0 I 

100 


I 10% 1 

1% 


1 $25,000 1 

$20,000 


1 5% I 

1% 


1 $62,500 1 

$50,000 


1 150 [ 

1.50 


I 0,020 [ 

0 002 


I 00 I 

00 


I 100% J 

100% 


] 100 I 

100 


I 25.0 I 

25 0 


I 100% j 

100% 


1 $500 1 

$500 


1 100% I 

100% 


I 10,0 I 

100 


I 3% I 

0% 


I $1,250 1 

$1,000 


1 1% 1 

0% 


1 $62,500 1 

$5,000 


1 150 I 

1 50 


I $100,000 1 

$100,000 


1 50 1 

$0 


1 ^ 1 

4 



0% 



0 00 




CEA Guru 
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title 

engine MODEL Fwww-»<-yyy 
TASK/ECP 


STAf^DARD HISTORY FILE 


(O) (P) 

No of Available Mod Months 















title 

ENGINE MODEL Fwww-xx-yyy 
TASK/ECP 


F-zz 


CEA VERSION 2 1 


12/9/94 
CEA Guru 



ASSUMPTIONS 

a) Incorporation in Production engines will begin in 

b) Number of engines produced with this change is 

c) Number of spare units incorporating this change is 

d) Modification of operational engines can begin in 

e) Incorporation of this change in operational 
engines will be accomplished by —> 

f) Total kits installed out of total 
engines not modified in production is 

g) Total engines lost to attrition is 

h) Total engines retired unmodified is 
1 } Estimated yearly flying hours 



May 

1993 



550 



0 


Aug 

1993 

1st Opportunity 

at 

Depot 

0 

of 

0 


50 

0 


240 EFH / Year 
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APPENDIX D: SIMAN SIMULATION PROGRAM 


The first page of this appendix contains the SIMAN Model 
Frame for the simulation program used to generate the Weibull 
data. The process involved in the program is full explained in 
Chapter IV. 

The second page contains the SIMAN Experimental Frame used 
in conjunction with the Model Frame to specify the specific 
parameters needed to run the program. The specifics of the 
experimental frame are fully explained in Chapter IV. 
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Begin,Yes,Yes; 

/Model for the Simulation of a Weibull Distribution relative 
;to the Aircraft Component Improvement Program 

; This is the Model portion of the SIMAN program. 

; Create Entity Module 

CREATE,10,0:,1; ' 

; Creates a distribution of failures of engine components 

; ten engines at time zero, max of one batch 

ASSIGN:X(1)=X(1)+1; 

ASSIGN:ENGNUM=X(1); 

; Assigns a sequential number 1 to 10 to each engine 

Q1 ASSIGN: Arrtime=TNOW; 

DELAY:WEIB(50,5); 

; Delays each engine for the operating period, end of delay 

; is the failure. 

QUEUE,RepairQ,• 

SEIZE:Repairman,• 

DELAY: 0; 

RELEASE: Repairman; 

; These steps represent the repair process, delay of zero to 

; simulate assumptions of CEAMOD program 

TALLY:ENGNUM,INTERVAL(Arrtime); 

; This step collects the data on the time of failure 

DELAY:0:NEXT(Q1); 

END; 
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BEGIN; 

PROJECT, Thesis, Glenn R. Cook; 


ATTRIBUTES: 

RESOURCES: 
QUEUES: 
TALLIES: 


This is the experimental frame of the SIMAN program 
it is designed to allow changes in the data for different 
program runs. 

Arrtime: 

ENGNUM; 

Repairman; 

RepairQ; 

1, Time for Engine 1, "Engl.dat": 

2, Time for Engine 2, "Eng2.dat": 

3, Time for Engine 3, "Eng3.dat": 

4, Time for Engine 4, "Eng4.dat": 

5, Time for Engine 5, "Eng5.dat": 

6, Time for Engine 6, "Eng6.dat": 

7, Time for Engine 7, "Eng7,dat": 

8, Time for Engine 8, "Eng8.dat": 

9, Time for Engine 9, "Eng9.dat": 

10, Time for Engine 10, "Engl0.dat"; 

REPLICATE, 10,0,240,N,Y,0; 

This indicates that the program will run once for 2400 
units of time. 
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appendix E; SIMAN OUTPUT REPORTS 


This appendix contains the output reports generated by the 
SIMAN program. The author has chosen to only include examples 
from the simulation runs, thus this appendix only has the reports 
from the runs with an eta value of 20. The reports are fully 
explained in Chapter IV. 
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SIMAN IV - License #9050352 
Naval Post-Graduate School 


EXPO (20) 


Project: Thesis 

Analyst: Glenn R.Cook 


Run execution date : 
Model revision date: 

11/27/1994 

11/27/1994 

Replication ended at time 

: 2400.0 




TALLY VARIABLES 


Identifier 

Average 

Variation 

Minimum 

Maximum 

Observations 

Time for Engine 1 

19.351 

.94969 

.20862 

83.182 

122 

Time for Engine 2 

17.520 

. 94349 

.35883 

73.424 

135 

Time for Engine 3 

19.741 

1.0464 

.12021 

103.66 

121 

Time for Engine 4 

21.066 

.95621 

.17511 

109.84 

113 

Time for Engine 5 

19.186 

.92286 

.09282 

96.620 

125 

Time for Engine 6 

19.296 

.83613 

.34302 

84.389 

124 

Time for Engine 7 

20.700 

.95878 

.27664 

116.35 

115 

Time for Engine 8 

20.483 

1.0106 

1.1141 

113.80 

116 

Time for Engine 9 

18.795 

.89282 

.10558 

90.237 

127 

Time for Engine 10 

21.799 

1.0095 

.07690 

108.14 

109 


Run Time: 0 min(s) 2 sec(s) 
Simulation run complete. 


Mean = 19.79 Hours 
Variance = 355.7 Hours 
Fail/Eng/Year = 12.07 
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SIMAN IV - License #9050352 
Naval Post-Graduate School 

Weibull(20,1.5) 

Project: Thesis Run execution date : 11/27/1994 

Analyst! Glenn R.Cook Model revision date: 11/27/1994 

Replication ended at time : 2400.0 

TALLY VARIABLES 

Identifier Average Variation Minimum Maximum Observations 


Time for Engine 1 17.710 
Time for Engine 2 17.271 
Time for Engine 3 16.077 
Time for Engine 4 19.058 
Time for Engine 5 16.727 
Time for Engine 6 19.110 
Time for Engine 7 19.280 
Time for Engine 8 17.493 
Time for Engine 9 19.460 
Time for Engine 10 17.619 

Run Time: 0 min(s) 3 sec(s) 
Simulation run complete. 


Mean = 17.98 Hours 
Varaince = 143.0 Hours 
Fail/Eng/Year = 13.31 


64641 

1.2111 

52.406 

135 

69992 

.55646 

62.322 

138 

65722 

.29150 

53.101 

149 

60095 

1.4608 

58.326 

125 

64935 

1.3782 

51.809 

143 

66742 

1.4340 

64.694 

125 

71869 

.85501 

63.743 

121 

67730 

.59326 

55.313 

137 

64598 

1.1069 

59.900 

122 

69236 

.84961 

69.638 

136 
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SIMAN IV - License #9050352 
Naval Post-Graduate School 


Weibull (20,2) 

Project: Thesis Run execution date : 11/27/1994 

Analyst: Glenn R.Cook Model revision date: 11/27/1994 

Replication ended at time : 2400.0 

TALLY VARIABLES 

Identifier Average Variation Minimum Maximum Observations 


Time for Engine 1 18.233 .52208 1.4531 50.979 131 

Time for Engine 2 17.019 .54696 1.4297 41.855 141 

Time for Engine 3 18.124 .48846 3.5033 39.402 132 

Time for Engine 4 17.884 .54376 1.2400 51.982 133 

Time for Engine 5 17.613 .52332 2.0428 46.506 135 

Time for Engine 6 16.724 .49560 2.3522 42.024 143 

Time for Engine 7 17.715 .48766 2.9432 45.533 135 

Time for Engine 8 18.128 .49125 2.7711 42.892 132 

Time for Engine 9 17.708 .54230 .83887 48.240 134 

Time for Engine 10 17.488 .50650 1.3625 45.068 137 

Run Time: 0 min(s) 2 sec(s) 

Simulation run complete. 


Mean = 17.66 Hours 
Variance = 82.7 Hours 
Fail/Eng/Year = 13.53 
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SIMAN IV - License #9050352 
Naval Post-Graduate School 


Weibull (20,3) 

Run execution date : 11/27/1994 

Model revision date: 11/27/1994 

Replication ended at time : 2400.0 


Project: Thesis 

Analyst: Glenn R.Cook 


TALLY VARIABLES 


Identifier 


Average Variation Minimum 


Maximum Observations 


Time for Engine 1 
Time for Engine 2 
Time for Engine 3 
Time for Engine 4 
Time for Engine 5 
Time for Engine 6 
Time for Engine 7 
Time for Engine 8 
Time for Engine 9 
Time for Engine 10 


18 

.885 

.35075 

17 

.435 

.37529 

17 

.560 

.33070 

17 , 

. 641 

.38849 

17 . 

,803 

.33332 

17 , 

.500 

.37213 

18 , 

. 031 

.37175 

18 . 

094 

.36543 

17 . 

681 

.33839 

17 . 

921 

.33384 


2.4143 33.810 
3.3360 35.705 
5.2358 33.048 
4.1221 37.807 
6.5280 33.660 
3.1330 37.320 
3.4824 32.318 
5.6019 35.971 
4.8010 29.547 
4.1352 35.286 


127 

137 

136 

136 

134 
136 
133 

132 

135 

133 


Run Time: 0 min(s) 3 sec(s) 
Simulation run complete. 


Mean = 17.85 Hours 
Varaince = 40.4 Hours 
Fail/Eng/Year = 13.39 
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SIMAN IV - License #9050352 
Naval Post-Graduate School 


Weibull (20,5) 

11/27/1994 
11/27/1994 

Replication ended at time : 2400.0 

TALLY VARIABLES 


Project: Thesis Run execution date : 

Analyst: Glenn R.Cook Model revision date: 


Identifier 

Average 

Variation 

Minimum 

Maximum 

Observations 

Time for Engine 1 

18.216 

.21669 

9.0715 

27.425 

131 

Time for Engine 2 

18.339 

.26445 

6.9617 

29.079 

130 

Time for Engine 3 

18.397 

.23844 

7.1913 

27.333 

130 

Time for Engine 4 

18.336 

.20049 

10.339 

26.597 

130 

Time for Engine 5 

18.592 

.20915 

6.5764 

27.137 

128 

Time for Engine 6 

18.693 

.21916 

8.6234 

28.126 

128 

Time for Engine 7 

18.137 

.22710 

6.8289 

27.794 

132 

Time for Engine 8 

18.587 

.21210 

9.3200 

27.573 

129 

Time for Engine 9 

18.436 

.20319 

7.7681 

27.680 

130 

Time for Engine 10 

18.332 

.22248 

7.7532 

28.317 

130 


Run Time: 0 min{s) 3 sec(s) 
Simulation run complete. 


Mean = 18.41 Hours 
Variance = 16.6 Hours 
Fail/Eng/Year = 12.98 
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APPENDIX Ft SIMAN "DAT" FILE OUTPUT 


This appendix contains the output of a SIMAN "Dat" file. 

The first coluran is the cumulative clock time from the first 
failure to the total run of 2400 hours. The second column 
represents the actual times to failure for engine 1. The use of 
this data is fully explained in Chapter IV. 
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Time for Engine 
.13679400E+03 
.16538160E+03 
.30977810E+03 
.43940080E+03 
.55595120E+03 
.62420570E+03 
.64869380E+03 
.66302080E+03 
.69013670E+03 
.82652700E+03 
.86262820E+03 
.96482040E+03 
.97945980E+03 
.10163370E+04 
.11069960E+04 
.12023170E+04 
.15486620E+04 
.16954530E+04 
.17580750E+04 
.19088880E+04 
.21682330E+04 
.22424920E+04 
.23284450E+04 
-.lOOOOOOOE+01 


1 

.13679400E+03 
.28587570E+02 
. 14439650E+03 
.12962270E+03 
.11655040E+03 
.68254580E+02 
.24488100E+02 
.14326970E+02 
.27115910E+02 
.13639030E+03 
.36101140E+02 
.10219230E+03 
.14639340E+02 
.36877080E+02 
.90659120E+02 
.95320680E+02 
.34634510E+03 
.14679170E+03 
.62621460E+02 
.15081300E+03 
.25934500E+03 
. 74259520E+02 
.85952880E+02 
.lOOOOOOOE+01 
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APPENDIX G; 

GRAPHS OF ANNUAL NUMBER OF FAILURES 











































EXPO (500) 
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3 
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4 
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5 

4 

6 

2 

7 

5 

8 

6 

9 

8 

10 

5 

More 
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WEIBULL (500,1.5) 



WEIBULL (500,2) 



Year 

Frequency 

1 

1 
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6 

3 
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6 
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10 

6 
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1 

0 
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4 

8 
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6 

4 
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7 
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More 
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I FREQUENCY FREQUENCY 


APPENDIX H; 

FREQUENCY OF FAILURE GRAPHS 


EXPO (20), FREQUENCY OF FAILURES 



0 2 4 6 8 10 12 1 4 16 18 20 22 I 

NUMBER OF FAILURES PER YEAR 


WEIBULL(20,1.5). FREQUENCY OF FAILURES 



0 2 4 6 8 10 12 14 16 18 20 22 

NUMBER OF FAILURES PER YEAR 


WEIBULL (20,2), FREQUENCY OF FAILURES 



NUMBER OF FAILURES PER YEAR 
WEIBULL (20,3), FREQUENCY OF FAILURES 




WEIBULL (20,5), FREQUENCY OF FAILURES 
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WEIBULL (20,1.25), FREQUENCY OF FAILURES 



0 2 4 


number of failures per year 




WEIBULL (20,1.75), FREQUENCY OF FWLURES 

t 


s 

QC 


10 



0 — 


Sill 



«= T o ^ 2 2 ? 8 


O <N 


number of failures per year 
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EXPO (50), FREQUENCY OF FAILURES 



2 3 4 5 6 

NUMBER OF FAILURES PER YEAR 


25 

20 

I” 

2 10 


WEIBULL 




i=j±J 



i 3 * 5 6 7 8 9 10 I 

number of failures per year 


WEIBULL (50.2), FREQUENCY OF FAILURES 



, 2 3 * 5 6 7 8 

NUMBER OF FAILURES PER YEAR 


9 10 


WEIBULL (SO.J), FREQUENCY OF FAILURES 



2 3 4 5 6 

number of failures PER YEAR 


WEIBULL (50.5), FREQUENCY OF FAILURES 
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EXPO (100), FREQUENCY OF FAILURES 



1 2 3 4 5 6 

NUMBER OF FAILURES PER YEAR 


30 



1 2 3 4 5 6 

NUMBER OF FAILURES PER YEAR 


WEIBULL (100,2). FREQUENCY OF FAILURES 



0 1 2 3 4 5 6 7 

NUMBER OF FAILURES PER YEAR 


WEIBULL (100.3). FREQUENCY OF FAILURES 



0 1 2 3 4 5 6 7 

NUMBER OF FAILURES PER YEAR 


WEIBULL (100.5), FREQUENCY OF FAILURES 



0 1 2 3 4 5 6 7 

NUMBER OF FAILURES PER YEAR 
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EXPO (500), FREQUENCY OF FAILURES 
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APPENDIX I: FINANCIAL SPREADSHEET OUTPUT 


This appendix contains examples of the output from the 
financial spreadsheet used to determine potential cost savings. 
The author has chosen to include only those with a comparison to 
the mean or eta value of 20 hours. The output from these 
spreadsheets is fully explained in Chapter V. 
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